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Abstract 

Using the improved quantization technique to the mini-superspace approximation 
of loop quantum gravity, we study the evolution of black holes supported by a 
cosmological constant. The addition of a cosmological constant allows for classical 
solutions with planar, cylindrical, toroidal and higher genus black holes. Here 
we study the quantum analog of these space-times. In all scenarios studied, 
the singularity present in the classical counter-part is avoided in the quantized 
version and is replaced by a bounce, and in the late evolution, a series of less 
severe bounces. Interestingly, although there are differences during the evolution 
between the various symmetries and topologies, the evolution on the other side of 
the bounce asymptotes to space-times of Nariai-type, with the exception of the 
planar black hole analyzed here, whose T-i?=constant subspaces seem to continue 
expanding in the long term evolution. For the other cases, Nariai-type universes 
are attractors in the quantum evolution, albeit with different parameters. We 
study here the quantum evolution of each symmetry in detail. 
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1 Introduction 



It has long been thought that any viable theory of quantum gravity should be able 
to resolve the singularities which afflict the classical theory. Of primary importance 
are the space-like singularities which are manifest in many cosmological models and 
black hole interiors. It is hoped that a quantum theory of gravity not only replaces 
the singularity with something more benign, but also allows us to answer the question 
of what happens beyond the corresponding singular geometry. 

Today, one promising theory of quantum gravity is loop quantum gravity (LQG). 
(For a review, the interested reader is referred to pQ- [5] and references therein.) Loop 
quantum gravity has already been successful in answering the question of the origin 
of black hole entropy [6]-|18]. In the paradigm of LQG this entropy has been shown 
to arise from counting the number of quantum micro-states which yield a macroscopic 
black hole of fixed area. By setting the Immirzi parameter to one specific value, the 
area/4 law is recovered for all black holes. The sub-leading contribution to the entropy 
has been found to be genus dependent [18], which is in agreement with studies done 
utilizing non-LQG techniques [21], [22], [23] . There has also been progress on this 
front from the spin- foam approach [19j . [20] . 

Regarding the issue of singularity resolution, there has been much work done in 
the framework of loop quantum cosmology (LQC) [M] and, more recently, with the 
Schwarzschild black hole [25]- [32]. In these studies, the full quantum theory is techni- 
cally formidable to deal with and therefore one resorts to first symmetry reducing the 
system before quantization, in the hopes that the essential features in the full quan- 
tum system remain intact in the reduced analysis. Studies so far, both in LQC and 
black hole interiors, indicate that the singularities are replaced by quantum bounces 
and therefore gravitation becomes repulsive at extremely high energies. 

In the black hole sector, detailed studies to date have been performed only in 
the case of the Schwarzschild black hole. As mentioned above, the Schwarzschild 
singularity is replaced by a bounce in the quantized version. Aside from the elimination 
of the singularity it is interesting to ask what happens beyond the bounce, as it is 
now possible to evolve the equations on the other side of the classical singularity. 
For the Schwarzschild black hole the evolution is not symmetric and one does not 
simply obtain a time-reversal of the original black hole interior on the other side of 
the bounce. It has been shown that instead, the solution asymptotes to a second 
non-singular black hole [30J in the case where the "length" of the holonomy path is 
taken to be constant (the "fixed-delta" scheme). In the scenario where the holonomy 
path is not constant (the "improved quantization scheme") the evolution tends to 
a Nariai space-time [30J. We will discuss the two schemes in more detail below. A 
particularly in-depth analysis of various evolution schemes, comparing their strengths 
and weaknesses is in [31] . 

We wish to study here other black holes and ask the question of whether or not 
the singularity is avoided and, if so, what is the subsequent evolution of the geometry 
beyond the classical singularity. We choose to study the A black holes. That is, 
black holes supported by a cosmological constant. The reason for this is that the 
admission of a cosmological constant allows for the existence of black holes of differing 
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symmetries and topologies in general relativity. For example, one may have planar, 
cylindrical, toroidal and higher genus black holes |21j . [23], |33j-|40|. (An interesting 
analysis of super-symmetric anti-de Sitter horizons may be found in [H].) Such a 
study is of interest for several reasons. The presence of a cosmological constant may 
be expected to play a non-trivial role on the large-scale evolution and therefore could 
affect the late-time evolution of the black hole interior (although it is not expected to 
play a significant role on the scale of the bounce itself). Secondly, if the singularities 
in these geometries are all avoided, then it will be evidence for the indication that 
space-like singularities are generically resolved within loop quantum gravity, and not 
just a product of the spherical symmetry or topology. Since the field of quantum 
gravity has a notorious lack of experimental guidance, it is the internal consistencies 
such as this which may ultimately decide if a theory is viable. Finally, the study of 
the final fate of the black holes is also interesting in its own right. In this paper we 
study the evolution and late-time behavior of A black holes of various geometries and 
topologies. We find that there are qualitative and quantitative differences during the 
evolution. However, in all cases the classical singularity is avoided and, at late time, 
the evolution tends to Nariai-type universes with the exception of the particular planar 
metric studied here. The values of the parameters in the asymptotic space-time and, 
of course, the symmetry of T = constant surfaces, depends on the initial space-time. 
Although the effective quantum theory is an approximation to the full theory, there is 
evidence in the case of Loop Quantum Cosmology that the effective theory similar to 
that presented here provides an accurate description of the quantum dynamics [32], 

This paper is organized as follows: In section 2 we review the formalism of the mini- 
superspace approximation for black holes in loop quantum gravity. In the sub-sections 
we begin by briefly reviewing the previously studied Schwarzschild case followed by 
detailed studies of the Schwarzschild-(anti)deSitter black hole, the planar black hole, 
and finally, the higher genus toroidal and hyperbolic black holes. Concluding remarks 
are provided in section 3. 



2 Quantum evolution 

In the canonical formulation of loop quantum gravity, one begins with a 3 + 1 decom- 
position of space-time and then re-writes the resulting action in terms of the Ashtekar 
variables; namely, a generalized SU(2) connection A l a and a densitized triad E°^, which 
is related to the three-space metric, q a b , via 

qq ab = I-.", F h / j . (1) 

These two variables then become the configuration and momentum variables respec- 
tively to be quantized, resulting in what are sometimes known as the quantum Einstein 

throughout this manuscript, units are chosen such that G=c=h — 1. 
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equations: 





d a A l b — d\,A l a + e l ■ k A 3 a A k b , and 7 the Immirzi 



5^E b jK a b (Kai, being the extrinsic cur- 



These equations are the quantum analogs of the now well known Gauss, vector and 
Hamiltonian (or scalar) constraints respectively. As mentioned in the introduction, 
the technical difficulty in finding and interpreting solutions to these equations are 
quite daunting. Therefore, in order to make some progress, one often first deals with 
a classical model, reducing the system to the desired symmetry, and then studies the 
quantum evolution of the reduced system. In the Schwarzschild interior and in LQC, 
the resulting Hamiltonian evolution leads to a difference equation for the quantum 
states, which is well behaved at the classical singularity (see, for examples, [29], [45]). 
However, enough progress has been made that there are now reasonably well defined 
prescriptions for generating classical Hamiltonians which are thought to encode the 
quantum dynamics of symmetry reduced systems [26]. We exploit this prescription 
here in the study of A black holes. The general scheme is as follows: One begins with 
a classical pair (A l a , E b k ) which a priori respect the desired spatial symmetry. One 
then partially gauge fixes the system by satisfying the Gauss constraint and writes the 
Hamiltonian constraint in terms of the remaining variables. In the resulting constraint 
the classical SU (2) connection components are replaced with effective holonomies of 
the connection (this is the effective quantum prescription) and the resulting quantum 
corrected Hamilton equations are solved for the connection and triad components. All 
the systems studied have the near horizon as their initial data surface. In this way we 
avoid large quantum effects at the horizon. This is similar to the procedure used in 
[30] as well as in [31J (more on this below). 

There have been several proposals for the effective quantum prescription to be used 
in such an effective Hamiltonian, mainly inspired by the much more studied systems 
of loop quantum cosmology. It has been shown that one may mimic the quantum 
dynamics by making the substitution aj — ► sin(a/<5) /5, where 5 represents the length of 
the holonomy path of the aj component of the connection in the principal (coordinate) 
directions. The early studies took 5 to be constant valued. However, it was later shown 
in LQC that using a fixed 5 prescription does not yield a good semi-classical limit [42J. 
If however the holonomy path- length 5 is made to depend on the momenta, Ej, a good 
semi-classical limit can be achieved. There are several candidates for this improved 
prescription and we use here the analog of the u Jl'" prescription of LQC [46J (details 
are provided below). We should mention before proceeding that one needs to use 
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caution when interpreting results within the mini-superspace approximation. There is 
no guarantee that this level of approximation always mimics the full quantum theory. 
An interesting recent analysis of this is in [47] . 



2.1 A brief review of the Schwarzschild interior 

Here we present a very brief overview of the previously studied Schwarzschild black 
hole. For full details the interested reader is referred to the pioneering papers [25] - 
|30j. One may choose a spherically symmetric connection and triad pair as 



A =o, 3 tr dR + {a\Tg + a 2 T^) d9 + (air^ — a 2 T$) sin 9 d<p + tr cos 6 dcj) , (5a) 
E =E 3 t r sin 6— + {E lTg + E 2 T (j> ) sin 6— + (E x t$ - E 2 t 9 ) ^- . (5b) 

Here the r represent the standard Pauli su{2) basis and aj and Ej are functions 
of the interior time only. The classical Gauss constraint yields only one non-trivial 
condition: 

2sin6>(ai£ 2 - E x a 2 ) = 0, 

which is customarily satisfied by setting a\ = E\ = 0. Setting E\ = still leaves the 
triad non-degenerate and it can be checked that the classical vector contraint is also 
satisfied. The classical Hamiltonian constraint now reads as: 

/N 
—^=[(a 2 2 +j 2 )E 2 + 2a 3 a 2 E 3 ]dRdn. (6) 

The R integral is performed over a finite length, Lq, and the resulting factor of L$ is 
absorbed into a re-scaling of a 3 and E 2 yielding: 

H = - N [{a 2 2 + 7 2 ) E 2 + 2a 3 a 2 E 3 ] . (7) 
27^ v h, 3 



In the lapse function, N, was taken as N = 7 and the Hamilton's equations 



v2 a 

(E 2 = {E 2 , H} = —dH/da 2 etc.) yield the following equations of motion: 



If 7 2 \ ■ If 7 2 

En — — — I CLo H ) , Eo — — I Eo — — kEo 

2 \ a 2 J 2\ a\ 

a 3 = — 2a 3 , E 3 = 2E 3 , 

with the dot representing the derivative with respect to the Hamiltonian affine pa- 
rameter, T. 

The solution to these equations yields, after a re-definition T — > T' := e T : 



»2=±7\/^-l, E 2 = Ef^/2MV -T<\ 



V 
T 



«3=T74 0) ^. E 3 = ±T' 2 , (8) 



5 



with and M constants. By comparing this solution with (pQ) it can easily be seen, 
after absorbing E% in a coordinate re-scaling, that the densitized triad is compatible 
with a metric which admits the following well-known line element: 

(2M S\ ,r,2 , r r /2./ ) 2 , rr>2 -2. 



da 2 = ^— - 1 J dJF + r W + T 2 sin 2 9 d<f) 2 . (9) 

This, of course, is the spatial line element corresponding to the Schwarzschild interior 
solution. 

The analogous quantum system has been studied in, for example, [29] where it was 
shown that a discrete eigenvalue flow equation arises. In a different approach [27], [30J, 
an effective Hamiltonian was utilized which encodes the quantum corrections. In both 
types of studies it has been demonstrated that the classical Schwarzschild singularity is 
avoided. In the effective Hamiltonian scheme, the quantization prescription described 
previously yields [30j : 



H 



2j sin(a 2 52 



sm 2 (a 2 d 2 ) , 2 \ p , sm(a3<5 3 ) sm(a 2 <5 2 ) 
75 rl )E 2 + Z E 3 



(10) 



where the <5j represent the length of the holonomy path in the I-th (angular or radial) 
direction. In the limit that the Sj's vanish, one recovers the classical Hamiltonian. 
The resulting equations of motion may be found in [27], [30j or in the A — ► limit of 
the following section (albeit in a different gauge). 

Studies where the fixed 5 prescription (5 2 = 5 3 =const.) has been used have found 
that the singularity is replaced by a quantum bounce, and that on the other side of 
the bounce the system evolves to another (non-singular) black hole |27| . [30] , In the 
non-fixed 5 prescription, where the <5's are functions of the triad (the scheme we utilize 
below), it was shown that the evolution on the other side of the bounce asymptotes 
to the Nariai space-time [30] . 

2.2 The Schwarzschild- (anti)deSitter interior 

The quantization of the Schwarzschild-(anti)deSitter interior is very similar to the 
above Schwarzschild system. We present it here mainly for completion and as a segue 
to the other A black holes studied below. 

The interior Schwarzschild-(anti)deSitter line element is given by: 



dT 2 /A„o 2M 

_ I 2M 
.3"* T 



ds2 = - TA^WTTj + I I T " + dR 2 + T 2 d9 2 + T 2 sin 2 9 d<p 2 , (11) 



where A is the cosmological constant and M the black hole mass. Due to the spherical 
symmetry, the form of the connection and triad are taken as in (|5ap and (|5b[) and the 
Hamiltonian constraint is therefore 

H = - [{a 2 + 7 2 ) E 2 + 2a 3 a 2 E 3 - AE 2 E 3 ] . (12) 

ZYVE 3 
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We choose a different gauge (choice of N) here, however, and pick the lapse function 

as 

7 2 v^ 



N = , (13) 

which amounts to working directly in Schwarzschild (T-domain) coordinates. That is, 
with this choice, and utilizing ([1]) and (jlip . the classical Schwarzschild-(anti)deSitter 
line element is given by 



ds l = ——-ydT + -^dR z + E 3 d9 z + E 3 sin z 9 d<f> z , (14) 
£•2 E 3 

(classically 7 can be set equal to unity without loss of generality) with: 



From these expressions it can be seen that the classical singularity occurs at T = 
where both £3 = and E2 = 0. The horizon(s) are present where only E2 = 0. One 
advantage of this choice of ./V is that, as we shall show below, the asymptotic form of 
the metric for large negative T is easier to find since E2 approaches a linear function 
of T. 

We now go over to the effective quantum system. In the current gauge the corre- 
sponding effective quantum Hamiltonian is given by 



H = - 1 - 
2 



sin 2 (a 2 <5 2 ) 2 \ , sin(o 3 <5 3 ) sin(a 2 <5 2 ) #3 A „ 

+ 7 J +2 ^ JT^*- m 



(15) 



We now address the issue of choosing the 5's. As shown in a cosmological setting in 
[42j [46j and for black hole interiors in [30] , a reasonable semi-classical limit might not 
be obtained utilizing a fixed 5. Instead, one can relate the length of the holonomy 
paths (the 5's) to the classical area. For example, in the R — 9 plane, the proper area 
spanned by the holonomy loop is given by: 

A/j e = E 2 5 2 S 3 = A , (16) 

whereas, in the 9 — (p plane a similar relation is given by: 

k e4> = 5jE 3 = A . (17) 

The final equality in (|16p and (|17p comes from setting the area equal to the small- 
est area gap predicted from loop quantum gravity, which we denote as A. These 
relationships yield: 

02 = \TT, 3 = = , (18) 

V ^3 -^2 

which are to be substituted in (|15p to complete the quantum prescription. The re- 
sulting equations of motion are: 
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1 



«3 



E 3 E 2 A 3 / 2 



-2y/ E 2 E 3 Asm [ — J sin I — 



/ asVAA . / a 3 VA^E 3 ' \ . ( a 2 VA \ ( agv^v 7 ^' 

A2AC0S — -=- sin — — a3£/3Asin — ;==- I cos 



E 3 f 

+ a2E 2 Asm | — ;==- | cos 



IE 3 

+ AE 2 A^/E 3 ~ = 0, 



E 2 



E 2 



(19) 



■ 2£ 3 3/2 . a 2 VA\ a 3 VAVE~ 3 \ 
E 3 — sin — ==- cos = , 



E 2 



(20) 



a 2 



a 3 E 3 3 / 2 . ( a 2 VA N 

/Eg , 



A^o 2 



sm 



cos 



a 3 



£"2 



0, 



(21) 



'E 3 cos 



^2 



sm — + sm — ==- 

[ e 2 [ 



0, 



(22) 



where the over-dot denotes differentiation with respect to T. These equations, along 
with H = are evolved using a surface within the horizon as the initial data surface 
(the inner horizon in the case of positive cosmological constant). Near this surface, it is 
expected that the classical Schwarzschild-(anti)deSitter solution is valid and therefore 
provides the initial conditions for the subsequent evolution. This should be a valid 
approximation as long as the black hole is sufficiently large as the curvature will be 
small here. However, the horizon itself can eventually pick up quantum corrections, 
as is hinted at by the fact that 5 3 becomes infinite there, and the <5's arise from the 
effective quantization. The study by Chiou [H] discusses in detail the effects of the 
quantum correction on the spherical, A = horizon where the black hole singularity is 
resolved and joins smoothly to the diffused horizon of the consecutive black hole. We 
study the long-time evolution in scenarios both with positive and negative cosmological 
constant. 

In figured] we show a sample evolution with positive cosmological constant. From 
figure QJ,) it can be seen that T = 0, where the classical £"3 = heralds the presence of 
a singularity, the corresponding quantum E 3 is small but non-vanishing. In fact, the 
evolution can be extended beyond T = where it can be seen that a series of bounces 
occur (see figures [lb) and c) for different domains.) These are analogous to the series 
of black holes discussed in [31] where the analysis revealed a fractal-like structure for 
the spherical black hole without cosmological constant. 

A set of sample results for the case of negative cosmological constant is shown in 
figure [2j 



S 




Figure 1: The quantum evolution for a Schwarzschild-deSitter black hole interior, a) The 
triad component E3 plotted in the vicinity of the classical singularity, b) The triad component 
E 3 plotted from T = -350 to T w -70. c) The triad component E c z up to T = -1 x 10 s . 
d) The triad combination E\jE 3 which is equal to the classical metric component gun- e) A 
close-up of the previous plot, near the classical singularity, f) A comparison of the classical 
(black) vs. quantum (red) evolution of E 3 for positive T. g) A comparison of the classical 
(black) vs. quantum (red) evolution of E2/E3 for positive T. h) The triad component E-2- 
The inner horizon is located at T « 10.4 and the following parameters were used: 7 = 0.274, 
M = 5, Aw 0.001. 

From the two sets of figures, it can be seen that although there are some small 
qualitative differences in the evolution of some of the components, the singularity is 
avoided in both cases. The component £"3 never vanishes. The late time behavior of 
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Figure 2: The quantum evolution for a Schwarzschild-anti deSitter black hole interior, a) The 
triad component E3 plotted in the vicinity of the classical singularity, b) The triad component 
E 3 plotted from T = -350 to T w -70. c) The triad component E c z up to T = -1 x 10 s . 
d) The triad combination E\jE 3 which is equal to the classical metric component gun- e) A 
close-up of the previous plot, near the classical singularity, f) A comparison of the classical 
(black) vs. quantum (red) evolution of E 3 for positive T. g) A comparison of the classical 
(black) vs. quantum (red) evolution of E2/E3 for positive T. h) The triad component E^. 
The horizon is located at T w 1.7 and the following parameters were used: 7 = 0.274, M = 1, 
A « -0.1. 

the two solutions is similar, approaching a Nariai-type universe similar to the scenario 
studied with A = [30]. For reference, a Nariai-type universe in the large T limit, in 
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the coordinate system chosen here can be written as 

ds 2 = ^-dT 2 + b -^dR 2 + a d9 2 + a sin 2 9 dcj) 2 , (23) 

OqT 2 a 

with ao and 60 constants. All orthonormal Riemann components are constants in this 
space-time. Furthermore, the radius of two-spheres is y/oo- The true Nariai universe 
(with all constant orthonormal Riemann components equal) is recovered in the limit 
&o = 1- From the graph of S3, we can see that the radius of the two-sphere is of the 
order of the Planck length. In the asymptotic regime, we assume that classical general 
relativity holds and that it is therefore justified to use the Einstein field equations, 
although caution must be used given that some the dimensions are near Planck size. 

2.3 The planar black hole interior 

In this section we examine the evolution of black holes with planar symmetry. We 
begin with a general plane-symmetric connection and densitized triad |49j : 

A = (A x x t x + A y x t y ) dx + {A Y y T Y + A x y T X ) dy + A\t z d(, (24a) 
E = [E\r x + E\r Y ) A + [E\r Y + E\r x ) |- + , (24b) 

with £ the space-like coordinate corresponding to the time-like coordinate in the ex- 
terior domain of the black plane. All coefficients are functions of the interior time, T, 
only. 

Before continuing, there is one subtlety that needs to be addressed for the subse- 
quent black hole scenarios which is not present for the spherical black hole. In the 
previous case, the fact that the scenario studied consisted of pure (Lambda) gravity 
with imposed spherical symmetry was sufficient to guarantee that the system under 
study was the Schwarzschild-(A)dS black hole, at least in the purely classical case. 
Such a guarantee is not present, even in the classical case, for other symmetries. 
That is, there is no Birkhoff 's theorem on which we can rely on. Admittedly, in the 
quantized evolution even in the spherical case we do not have a uniqueness theorem. 
However, the fact that Birkhoff 's theorem holds (at least approximately) up until the 
quantum effects become strong should help insure that we are studying the analog of 
the Schwarzschild-(A)dS until the late-time quantum evolution takes over. To address 
this issue in the non-spherical cases, we test the classical evolution equations by com- 
paring their results in the classical case with the known exact solutions. We set a tight 
tolerance in the evolution equations so that the numerically evolved classical equations 
agree with the exact solution within the tolerance throughout the entire evolution. We 
then "quantize" the equations of motion and evolve the quantum system using this 
tolerance. 

The simplest planar black hole interior metric is given by [21] 

ds 2 = + (f + ^) d( 2 + T 2 (dx 2 + dy 2 ) , (25) 
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where k is a constant related to the mass-per-unit-area (recall A < here). The 
classical time coordinate T possesses the range < T < ?/ ^jr where the upper limit 



corresponds to the event horizon. This metric will be used as a consistency check for 
our evolution scheme as well as to provide the initial conditions for both the classical 
and quantum evolutions at the horizon. However, it should be emphasized that the 
momentum components of E and the configuration components of A are treated as 
free parameters spanning the minisuperspace of the system and the geometry is not a 
priori fixed, save for the initial data surface. 

The Gauss constraint, utilizing (|24ap and (|24b|) yields only the following one con- 
dition: 

- A Y X E\ + A X X E\ - A\E\ + A X E\ = . (26) 

Satisfying this constraint will also automatically satisfy the vector constraint. We can 
gauge fix the system in a way to satisfy this constraint and be compatible with the 
metric (125p (via ([I])) with the prescription: 

A x x = A Y y , E\ = E y Y , A Y X = 0, A x y = 0, E\ = 0, E\ = 0. (27) 

The equality of the two diagonal components is inspired by the planar symmetry of the 
system whereas setting the off-diagonal terms to zero satisfies the Gauss constraint!! 
(126h . Before continuing, we mention that although A x = A Y y and E x x = E y Y , we treat 
them here as independent degrees of freedom and therefore they are evolved as such. 
We will show below that with the appropriate initial conditions of A x (Tq) = A Y y (To) 
and E x x (Tq) = E v y (Tq), the equations of motion will preserve this equality for all 
subsequent time. 

With the above ansatz, the resulting classical Hamiltonian reads: 



H 



2N 



E\E\E^ Z 



E\E\A X A Y + E\E^A X X A\ + E? Y rf z A Y y A z c 



+ ^E* xE y Y E< Z: (28) 

where quantities have been re-scaled in order to absorb the length factors that arise 
from the truncated spatial integrations. We choose to work in a coordinate system 



compatible with ([25]) and therefore take^l N = ^ z . Utilizing ([T]) and (j25j) . this 
allows us to write the line element as 



M, o E x 2 



ds 2 = - -^dT 2 + -^dQ 2 + E C Z (dx 2 + dy 2 ) . (29) 
E*X E\ 



2 For the planar case the fiducial spin connection r ! a vanishes and therefore there is no inconsistency 
in choosing A % a as diagonal. 

3 An ambiguity exists regarding whether to take E x x or E y Y (or some combination) in this ex- 
pression. In this classical case, it makes no difference. We discuss the issue below in the quantum 
case. 
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This choice leads to equations of motion: 
2 



E x x 3 / 2 JE y Y 1 



E\A x x A Y y + E%A X X A Z 



+ 



E x x 3 / 2 y E y Y 
AE C z E y Y 



E X X V\ E V y 



E x x E y Y A x x A Y + E\E^A X X A\ + E\E^A Y y A\ 



0, 



(30) 



AYy + E x3 I 2 Je\ 



x 
1 



Y 



E\A x A Y y + E^AlA z r 



E X X 3 I 2 E\ 3 I 2 



E\E\A x x A Y y + E\E%A X X A\ + E y Y E%A Y A z 



AE^ 



E X Ey 



o, 



(31) 



F, x — 



e x x 3 / 2 Je v y 



E\E\A Y y - E X X E< Z A Z 



(32) 



E y 



Y 



E x x 3 / 2 \/ E y Y 



E\E\A X X + E y Y E< z A z 



(33) 



A z + 



E X x 3/ 2 JE y y 



[E\A X X A\ + E\A Y y A 



AJE y 
) A z^ V ■ 



x 
X 



(34) 



E' 



E x x 3 / 2 



E\ 



E x x EiA x + E y v EiA Y , 



(35) 



Although complex, these equations can actually be solved analytically. However, 
as mentioned above, we wish to test the computational evolution scheme and therefore 
we solve the above numerically, using the exact classical solution to set the values of 
the triad on the initial data surface just inside the horizon. The initial connection 
components are set from first solving the classical Hamiltonian equations for A z ^, as 

we know the classical E x x and E z from the form of the metric (|25p . In figure [3] we 
show the ratio of the numerical results to the exact analytic value and can see that 
the ratio is very close to unity throughout the evolution. It can be seen that near 
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a) 



b) 



2 4 



d) 



e) 



f) 



Figure 3: A comparison of the numerically evolved classical solution to the exact classical 
solution (denoted with subscripts "e"). We show here the ratios of a) E^ z /E^ ze , b) E x x / E x xe , 
c)A z ( /A z ce , d) A x x /A x xe . The ratio is very close to unity, indicating that the numerical solution 
is equal to the exact solution throughout the entire evolution. In e) and f) we show a close-up 
of the ratio of E V Y /E x xe and A Y y jA x xe respectively. These plots indicate that the equality 
of these components is preserved throughout the evolution, compatible with the ansatz (|27p . 
The horizon is located at T ps 12.5. 



the singularity as well, the ratio is acceptably close to unity. We are therefore indeed 
evolving the planar black hole metric of (|25p . It can also be seen from these figures 
that the equalities E x x = E y Y and = A Y y are preserved throughout the evolution. 

Next we study the corresponding quantum system. The quantization of the Hamil- 
tonian proceeds as in the previous spherical case. However, there are some subtleties 
which must be addressed first. The classical three-space metric is given by (|25p . 
However, when this metric is written in terms of the Ashtekar variable E, there is 
an ambiguity due to the equality of E x x and E y y . For example, all of the following 
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combinations are equivalent: 



ds 2 
ds 2 
ds 2 
ds 2 
ds 2 



-^ 2 dT 2 + ^- dC 2 + Efi z (dx 2 + dy 2 ) 



< 



j, » j^y 2 

~ dT2 + i?r dC ? + E ^ ( dx2 + d y 2 ) ■ 



2 1. -/ ., ' . v 

A dT 2 + 



+ E V y 2 



X + c Y d( 2 + £ C Z (dx 2 + dy 2 ) . 
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(E x x Ey 



A rJT 2 A 



ExEy ^ d( 2 + £ ? 7 (cix 2 + <iy 2 ) , 

j-,C, v / 



(36a) 
(36b) 
(36c) 
(36d) 
(36e) 



or any other combination which yields the same coordinate dependence for the metric. 
In the classical evolution, this ambiguity makes no difference to the subsequent evo- 
lution. However, in the quantum evolution there is a difference as the metric is used 
to define the length of the holonomy paths via the planar analog of (|16p and (|17p . It 
seems intuitive that one should pick one of the symmetric choices, similar to (|36c|) . 
(|36dp or (|36ep . In this case, it turns out that all of these symmetric schemes yield the 
same results and therefore the bulk of our results here utilize (|36ep . which yields the 
most simple Hamiltonian. The evolution utilizing (|36ap or (|36bp is, however, quali- 
tatively different from the symmetric schemes. Below we show a sample evolution of 
this scenario as well. The important point is that all scenarios tested removed the 
classical singularity but there are sometimes significant differences in the evolution. 
The quantization procedure is similar to the spherical case. Briefly, we define Ss 

by: 

(37) 



"»xy 



E C Z S 2 



A, A 



E- X E V Y 5 ( 5 X 



A. 



where, as before, A is the minimum area predicted by loop quantum gravity. This 
yields: 



A 



y/A\/E C 



E X Ey 



(38) 



The quantum Hamiltonian reads: 



H 



2N 



^'E^E^E^ 



E x E y sm(A x x S x )sin(A Y y S x 

X ' Sr. S r . 



c sin(A x x S x ) S in(A z ( S ( ) 
+ ^X^z J F 



f sm(A x x S x 
+ E Y E z 7 



sin (.4^ 6{) 



2 7 2 



+ ^tjE^E^Ei 



(39) 
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We choose the gauge 



N= r^— 1 , (40) 

jpx jpV 

E X E Y 



which corresponds to using the coordinate system in (|25p in the particular "momentum 
gauge" of (|36e|) . This leads to the following "quantum" equations of motion: 



ax E z 

x A 3 / 2 E y Y 3 / 2 E x v 5 / 2 



' x 



E C z E y Y A z c AE x x ^ 2 sin | ^ r 



E x x 2 E y Y VAsm I ^4 



sin 



A^AE^ 



'E X X Ey 



+ E y Y 3 / 2 yj E C Z E X X A Z C A sin 



A^VA 



cos 



A\\j AE^z 



E\E\ 



+ E y Y 2 E x x VAsm 



A Y y VA 



sin 



E x Ey 



COS 



A z JAEi 



E X Ey 



(41) 



E a 



'E\cos 



E y VA 



E y Y sin 




E x x E y Y sin 




(42) 
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2JE^ Z E\E\ l 



+ 2£ a x £ y y Asin 



cos 



+ aa 3 / 2 Je c z e x x e v y 



E* X E\ 



2E X X E C Z A Z C A sin 



cos 



a z Va^e ( z 



E\E\ 



4E y v 3/2 \/ E x x E ( z VAsm 



AlVA 



sin 



a z ^Je c z 



E* X E\ 



+ 2E y Y 3 / 2 ^/E*~A Y y Acos 



sin 



a z ( Va^e c z 



E X Ey 



2E y Y E ( z A z (; A sin 



cos 



a z c Va^e c z 



E X Ey 



(43) 



E x x E y Y VA 



E\ 



A*y/A 



x 



sm 



+ £^ y sin 



= 0, 
(44) 



where the over-dot denotes differentiation with respect to the parameter T. The y — Y 
equations are related to the x — X equations via the interchange of x — > y and X — ► Y 
and vice-versa. 

As before, we utilize classical values at a point near the horizon as the initial 
conditions. Results for one set of parameters is shown in figure IH 

Note from (|36cp that the classical singularity at T = corresponds to E z = 0. 
However, from figures E^, b, and c it can be seen that the evolution now proceeds 
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Figure 4: The quantum evolution for a planar black hole interior, a) The triad component 
E z plotted in the vicinity of the classical singularity, b) The triad component E z plotted 
in the vicinity of the bounce, c) The triad component Eg up to T = — 1 x 10 s . d) The 
triad combination (E x x E y Y ) / E^ z which is equal to the classical metric component g^. e) A 
close-up of the previous plot, near the classical singularity, f) A comparison of the classical 
(black) vs. quantum (red) evolution of E z for positive T. g) A comparison of the classical 
(black) vs. quantum (red) evolution of (E X X E V Y ) j 'E% for positive T. h) The triad component 
E x x . i) The connection component ratio /A* . j) The triad component ratio E Y /E X . The 
horizon is located at T = 1 and the following parameters were used: 7 = 0.274, n — 0.00075, 
A w -0.0022. 
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through T = without E z vanishing. In fact, E z does not vanish for any value 
of T but instead oscillates and then smoothly starts growing again as we go to large 
negative times. Therefore the singularity present in the classical theory is avoided 
in the corresponding quantum evolution. In figure [3H we plot the densitized triad 
combination E x x E y Y /E^ z , which corresponds to the classical metric component g^. 
From figure Hfe, which is a close-up near T = 0, it can be seen that this is also finite 
near the classical singularity. Figures Hf and g compare the quantum evolution of 
E^ z and E x x E y Y /E^ z = to the exact (classical) solution from T ~ to the near- 
horizon. It can be seen that the quantum quantities are almost coincident with the 
classical quantities in this regime as the two plots on graphs overlap significantly. The 
quantum effects therefore start to become important only very close to T = 0. Finally, 
figure Hh shows the evolution of the triad component E x x (= E y Y ) and figures i and j 
show the consistency of E Y = E x and = A% throughout the evolution. 

From Fig. 4c), E^ z appears to grow without bound as T becomes more negative. 
If this is the case, then the space cannot be Nariai-type like the other cases. 

In figure [5] we also show a sample evolution based on ([36a]) for comparison. It may 
be seen that again the singularity is avoided. However, with the same parameters as 
the previous evolution scheme, the rate at which the asymptotic solution is achieved 
is different than with the symmetric scenario. 




Figure 5: The quantum evolution for a planar black hole interior using a non- symmetric evo- 
lution scheme, a) The triad component E^ z plotted in the vicinity of the classical singularity, 
b) The triad component E z plotted for large negative T. c) The triad combination E x x 2 /E z 
which is equal to the classical metric component g^^ in the vicinity of the classical singularity, 
d) The triad combination E x x 2 / 'E z for large negative values of T. The horizon is located at 
T = 1. 

In closing this section we note that the results for the flat cylindrical black hole 
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should be similar to the planar case. The flat cylinder arises from an identification of 
one of the coordinates, x or y in the metric (|25p . 



2.4 The higher genus black hole interiors 

In this section we discuss the interior of higher genus black holes. A reasonably general 
line element describing such interiors is provided by [35] 

ds 2 = ~ 2_M_ d p + A T2 + { 2 -T-( 3 + | r2 ) dR2+T2 de 2 +T 2 csinh\Vd9) # 2 , (45) 

with A < 0. The two cases of interest here are the following: 

i) j3 = 0, d = 0, lim c = -j: In this case constant (T, R) surfaces are tori. 

d — >0 

ii) /3 = — 1, c = d = 1: In this case constant (T, R) surfaces are surfaces of constant 
negative curvature of genus g > 1, depending on the identifications chosen. 

An event horizon exists when (tjt- — /3 + f^ 2 ) = 0. Black holes with higher genus 
topology have also been studied in [21] and [34] . 

The connection and densitized triad, for the SO (2, 1) symmetry of the constant 
(T, R) subspaces, are found to be 

A- =o-2,tr dR + (a\TQ + a 2 T (t) ) d6 + (a 2 r e - a\r^) \fcsmh{\fd6) dcp 

+ T fiv / cv / dcosh(v / d6>)#, (46a) 

E = - E 3 T R ^smh(Vd6)— - {Erre + E 2 r^) ^smh(Vd6) — 

+ {E X Ti - E 2 T d ) A . (46b) 

In the limit d = 0, and lim c = jjj this also yields the symmetry required for the 

toroidal case. That such a pair is indeed invariant under the appropriate symmetry 
group can be confirmed via the existence of a function W such that 

C^A = DW, C^E = [E, W] , (47) 

with the operator D being the gauge covariant derivative, d a W + [W, A a ] ( [ , ] is 
the su(2) Lie bracket). is the Lie derivative in the direction of the Killing vectors 
generating the symmetry, which for SO(2, 1) are: 

£ x := sm(<f>)dg + cos(0) cot h (0)9,/, , £ 2 := — cos(4>)dg + sin(</>) coth(#)<90 , £ 3 := . 

(48) 

It can be verified that {W = cos(</>)/ sinh(0)rij and 2 W = sin(0)/ sinh(0)TR satisfy 
(|47p . In the planar case, the situation was trivial, as setting W = satisfied (|47p 
which is then simply a statement that ^4 and E were independent of x and y. 
The Gauss constraint again yields a single condition: 

2 v / csinh fVdff) [a 2 £i - aiEfe] = , (49) 
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which will be satisfied by setting a± = E\ = 0. The vector constraint will also be 
satisfied by this choice. The resulting classical Hamiltonian is given by 



H 



N 



cos 



h(Vd6»i) 



[(a 2 2 - 7 2 d) E 2 + 2E 3 a 3 a 2 - AE 3 E 2 ] , 



(50) 



where 9\ represents the upper limit of the integration. In order to deal directly in 
the coordinate system of ([4"5"|). we choose the following: 



N 



''VdVE~ 3 



yfc cosh(\/d^i) — 1 



E 2 



(51) 



At this stage, the quantum prescription is similar to the scenarios in the previous 
sections. That is, the length of the holonomy paths are taken as 




E 2 



which is appropriate as (|45p admits proper areas of 

*rs = E 2 5 2 5 3 = A , A e(f> = 5 2 2 E 3 



The above Hamiltonian is then quantized via 



« 2 



sin(a 2 52 



03 



sin(a 3 ^ 3 
63 



(52) 



(53) 



(54) 



At this stage we should note that the particular choice of N given in (|5ip possesses 
an interesting advantage over other choices of the lapse function. Note that with this 
choice of N, the term in the Hamiltonian which depends on the parameter d decouples 
from the degrees of freedom. That is, with this choice, the evolution equations for 
the S-adS, toroidal and higher genus scenarios are similar. The difference between 
these cases then lies only in the initial conditions, as, for example, both the initial 
conditions for E 2 and a 2 involve the parameter /3, which is different for the three cases. 
We therefore do not present the equations of motion again here but simply show the 
results of the evolutions for both the toroidal and higher genus cases in figures [6] and 
[7] below. 

The higher genus hyperbolic evolution was the most computationally intense and 
therefore the evolution was halted at T = —10000 for the g > 1 case. This evolution is 
also the "noisiest" and produced the highest frequency of the resulting oscillation, as 
can be seen from figure [71 If these oscillations are a genuine product of the evolution, 
then it is perhaps more appropriate to state that the higher genus case approaches an 
asymptotic state which oscillates about a Nariai space-time, as the average of E 3 is 
seemingly approaching a constant. 

The late-time evolution of the torus was also a Nariai- type universe, similar to (|23p . 
with the appropriate toroidally compatible 2-space symmetry. We performed several 
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Figure 6: The quantum evolution for a toroidal black hole interior, a) The triad component 
£3 plotted in the vicinity of the classical singularity, b) The triad component E3 plotted from 
T = — 2 to T « 0. c) The triad component E 3 up to T = — 1 x 10 8 . d) The triad combination 
E2/E3 which is equal to the classical metric component gnu- e) A close-up of the previous 
plot, near the classical singularity, f) A comparison of the classical (black) vs. quantum (red) 
evolution of E 3 for positive T. g) A comparison of the classical (black) vs. quantum (red) 
evolution of E\jE 3 for positive T. h) The triad component E^- The inner horizon is located 
at T w 3.8 and the following parameters were used: 7 = 0.274, M = 1,A» —0.1. 

toroidal evolutions and noticed that in all of them, the frequency of the late-time 
oscilations was also greater than in their spherical and planar counter-parts. 

In [50] Bohmer and Vandersloot show that the Schwartzschild interior solution is 
a stable solution of the LQG equations. We suspect that all of the interior solutions 
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Figure 7: The quantum evolution for a higher-genus hyperbolic black hole interior, a) The 
triad component E 3 plotted in the vicinity of the classical singularity, b) The triad component 
E 3 plotted from T = -255.5 to T w -254. c) The triad component E 3 up to T = -10000. 
d) The triad combination E\jE 3 which is equal to the classical metric component gun- e) A 
close-up of the previous plot, near the classical singularity, f) A comparison of the classical 
(black) vs. quantum (red) evolution of E 3 for positive T. g) A comparison of the classical 
(black) vs. quantum (red) evolution of E2/E3 for positive T. h) The triad component E-2- 
The inner horizon is located at T « 6.2 and the following parameters were used: 7 = 0.274, 
M= 1, A « -0.1. 

we have discussed will likewise be stable. 
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3 Concluding remarks 



In this manuscript we have studied the evolution of the Schwarzschild-deSitter, Schwarz- 
schild-anti deSitter, planar, toroidal, and higher genus hyperbolic black holes in the 
mini-superspace approximation to loop quantum gravity. We have utilized the im- 
proved quantization scheme where the length of the holonomy paths are not constant, 
but instead are related to the proper area dicated by the corresponding classical den- 
sitized triad. We have found that in all cases, the classical singularity is avoided and, 
at some point during the evolution a "quantum bounce" occurs, usually followed by a 
series of less severe bounces. The exact details of the bounces and evolutions depend 
on the symmetry but the indication from all of these evolutions is that Nariai-type 
solutions seem to be generic attractors for the asymptotic space-time on the other side 
of the bounces for most cases studied. This is also an indication that loop quantum 
gravity is consistent in replacing classically singular space-times with non-singular 
ones, regardless of symmetry, at least at the level of approximation employed here. 
It would be interesting to improve on this scheme even further to see which of these 
features remain intact as one lessens the degree of approximation. Progress along this 
front is already being made. For example, see discussions in [47] and [51j . 
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